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Original English Summary. - A systematic method of constructing potentials, for which 
the one-variable Schroedinger equation can be solved in terms of the hypergeometric (HGM) 
function, is presented. All the potentials, obtained by energy-independent transformations of 
the HGM equation, are determined together with eigenvalues and eigenfunctions. A class of 
potentials derived from the confluent HGM equation is found by means of a limit process. 

To study theoretically the rotational and rotational- vibrational spectra of diatomic molecules, 
one often uses one-dimensional model potentials for which the solution can be expressed in terms 
of the HGM functions Q . The problem of the changes of variable in the HGM equation leading to 
the one dimensional Schroedinger equation has been studied several times in the past |l|, ^, [J, 
but the quoted authors have considered only the case of potentials given in explicit form Q. 
However, knowing the energy spectrum and the wavefunctions of the equation may prove useful 
for many applied problems whether or not the potential is given in explicit or implicit form. 

The general form of the transformations leading from the HGM equation to the Schroedinger 
equation is determined by the requirement (z = dz/dx) 

J\*t/~\ , 1 r, „i _ ;„2 



[z yi{z) + -{z,x} = Ar — 2MU(x) 



where 



_ (l-Ag)(l-z) + (l-Af)z + ( M 2 -lMl-z) 
I[Z) ~ 4z 2 (l-z) 2 ' {Z) 



lr The existence of potentials U(x), not explicitly depending on x, has been first posed in M (see also ff). 



and {z,x} is the Schwartzian derivative of z(x) with respect to x 



{z,x} 



z 6 z 
z 2 z 



(3) 



U (x) is the potential function; M is the reduced mass; E = J^j is the energy. When condition 
(1) is fulfilled the wavefunction is related to the HGM function as follows 

V[z{x)\ = (z'y^z^ (1 - z) h r 1 F(a, 0, r , z) . 



(4) 



Here 



Ao = 7- 1, 

Ai = a + P - 7, 

H = ft — a. 



(5) 



We shall assume that z(x) and therefore also {z, x} do not depend on the energy E. In this 
case, comparing the lhs and rhs of Eq.(l) one concludes that the parameters ^u 2 , Aq and \\ are 
linear in k 2 



f l_ M 2 = afc 2_ / 

) 1 — X 2 = c p k 2 -h p p = 0, 1 
and therefore z(x) fulfills the differential equation 

{z'fR(z) 



(6) 



Az 2 {l-zf 



where 



R(z) = a{z — l)z + cq(1 — z) + ciz 



(7) 



(8a) 



or in a slightly different form 



R(z) = az 2 + boz + cq = a(z — l) 2 + b\(z — 1) + ci 



(86) 



Introducing in the logarithmic derivatives of z' and z" that can be found using Eq.(7), 

one next gets the potential from Eq.(l) by means of Eq.(6) 



2MU[z(x)] = /^-D + Ml-^) + M + l + 

R 



a + 



a + (ci - Cq){2z - 1) 5 A 



z(z - 1) 



4R 



z 2 {l-zf 



R 2 



(9) 



where the discriminant A = b 2 — 4ac p = (a — c\ — cq) 2 — 4ciCq. 
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Thus, one will get a semiparametric (i.e., including the constant of integration of the differ- 
ential equation) family of potential curves, for one can always choose in the role of the three 
parameters, the scale factor, the origin of the coordinate x and the origin of the energy E. 

For the Schroedinger solution given by Eq.(4) the variable z is considered within the interval 
[0,1]. The corresponding interval for the variable x will be studied shortly. 

Since the case of the potential well of infinite depth is not of much interest, we shall ask the 
three-term quadratic polynomial R(z) to have no zeros for < z < 1, that is we shall consider 

R(z) > . (10) 

It follows from Eq.(7) that in the unit interval the function z(x) is monotone so that the corre- 
sponding transformation is single valued. Eq.(10) implies ^ 

c p > (11a) 



Integrating Eq.(7) gives 
a) A / 



±2x = 

p=0 



ijcp In | bp 2 y/ Cptp | bp J ^ ^ 



(12a) 



where t p = ^ z _^ ' ; 



b) A = 

l 

±2(x-x ) = J2(- 1 7V^ ] n\ z -P\ ■ ( 12b ) 
The transformation z(x) can be obtained in explicit form only for a few particular values of 
the parameters a and c p . One should notice the case when the zeros of R coincide with the 
singularities of the HGM equation. The resulting potentials have been already considered by 
various authors as follows 

a) Poeschl- Teller potentials M (i? = boz(l — z) or R = bo): 



h +l h!-< 3 



2Mb U(x) 



f + l sin 2 (x/y^) cos 2 (x/V^) ( 13 °' ) 



For A > 0, a > the coefficients bo and b\ should be of the same sign, i.e., 6q6i = (ci — Co) — a > 0. 



b) Rosen-Morse potentials || (R = cq) 



on* TTf \ h + In +2 ho-hx 1 / 

2Mc U(x) = + -^-thCx/v^) - l db 2 {x/ ^ (14) 



c) Manning-Rosen potentials QJ (R 



az 2 ) 



2MaU( X ) = + l^ctHx/^) ~ \^JT^ • ( 1 « > 

I 2 4 sn (x/y/a) 



The potentials (13-15) have the following general properties: the isotopic shift can be described 
by means of some changes of the parameters /, Hq and h±. This property does not extend to the 
whole family of potential curves given by Eq.(9). The scale transformation of the coordinate x 
leads in this case to an equation with five singular points, 0, 1, Zi, Z2, and oo, where z\ and Z2 
are the zeros of R. 

Let us find now the interval where x is defined for A / 0. For that, we first notice that the 
integral entering Eq. (12a) has no singularities, because from t p = yfa or 



U— + -, ^ N9 = y/a + 



z — p (z — p) 2 z — p 

A = would follow. This is why, performing z — > p (t p — > b p /2^/c^) in Eq. (12), we get the 
limits ±oo in x for c p ^ 0; if however one of the coefficients c p is nought then the corresponding 
choice of the coordinate origin for x as well as of the sign in Eq.(12a) can be done such that 

< x < oo . 

Finally, the variable x has a finite definition interval whenever both cq and c\ are zero. 

The case when only one of the c p coefficients, for example cq, is zero, is of special interest for 
the theory of two atom molecules, because for x — >• (z — » 0) the function 2MboU(x) goes to 
infinity as (/io + f ) z ~ 1 - F° r this potential, contrary to the Morse case, the Schroedinger equation 
is solved with regular boundary conditions for x = (see jH). 

If ho = | , the potential U(x) is finite at the origin and for x > (0 < z < 1) reads 

o/i/f rrr I M f(z Z l l ± hl + 1 , i , 3 (3*i - 1)(1 - *i) 5 ^i(l-^i) 2 
2Mfl[/[z(x)] = ^ + 1 + 4 (z - ^ " 4 (i - , l} 3 ( 6) 
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For x < it is natural to a priori determine the potential Eq.(16) to be symmetric U(x) = U (— x). 
Notice that for 

the potential Eq.(16) has two symmetric minima separated by a small barrier. 

Since the HGM solution which is irregular at z = 1 has a (1 — z)~ Xl (X\ > 0) behavior at 
that point, the integral 

/ m 2 dx = - J Rz A °- 1 (l-z) Al ~ 1 [F(a,/3,7;z)] 2 dz (17) 
Jo 4 jo 

is finite only if a (or (5) is equal either to a negative integer or zero: a = —n. The eigenfunctions 
are determined by the conditions 



dx 



n 'x=o = (18a) 



and 

*-(0)=0, (186) 
leading to a finite integral Eq. (17). The energy spectrum is determined by 

P + \ + yjh + 1 - Cl kj = y/f + 1 - akj , (19) 

where p = 2n for even levels and p = 2n + 1 for the odd ones, respectively. 

Let us clarify now under what conditions the potential Eq.(9) leads to a discrete spectrum. 
We shall consider c\ ^ 0. For a discrete spectrum to occur, the function (z')~Hfr should go to oo 
when z — > at least as z~z , if one takes ^ as the general solution of the Schroedinger equation. 
This condition is obviously true for cq ^ 0. If cq = the existence of the discrete spectrum is 
possible only for positive values of ho. The energy levels come from the equation 



2n 



+ 1 = y/f + 1 - akl - y^o + 1 - cokl -yjfn + l- Cl kl . (20) 



It is easy to see that the spectrum is bounded from above, if and only if both Co and c\ are 
not zero, and therefore, it follows from Eq.(20) that the discrete part of the spectrum has only 
a finite set of energy levels. Only the potential curve Eq.(13a), for which cq = c\ = 0, has an 
infinity of discrete levels. 



The eigenfunctions of the discrete spectrum are of the form 



f n = B n z^{\ - z) h 2 1 (Rz')^P^' X2 \2z - 1 



(21) 



HerePi Al ' Ao) (2z-l) are the Jacobi polynomials, whereas B n is a normalization constant given 



by 



B,, 



ci c a\ T(A + n + l)r(Ai +n + 1) 



+ 



(22) 



,Ai Ao nJ n\T(fi — n) 

(see the integrals 7.391 (1) and (5) in ||). In particular, the eigenfunctions VPn , corresponding 
to the potential Eq.(16), are obtained from Eq. (21) for c$ = 0, Ao = 

We notice that the potential curve Eq.(13a) displays one more important property: the Jacobi 
polynomials for neighbour energy levels are connected by simple recurrence relationships (the 
parameters Ao and Ai does not depend in this case on energy), and the matrix elements {m\z\n) 
is different from zero for m = n, n ± 1. 

To pass to the limit of the confluent HGM equation we make the following notations 



and 



°2 0-2,0-1, /, 01- 

a = -o, ci = -r H hc (6 = — 

T z T z T T 



f 



92 , 92 . 91 . , 
— oi ft i--ni r no 



z = tC ■ 



(23a) 
(236) 

(24) 



For t — ► we have 



where 



HO 



(CM = 1 

4(2 

^2 <?! 1-Ag 

4 4C 4C 2 



<*1 = hm r _*o[r(Ai - ^ 2 )] = g 1 - axk 2 . 
5 2 = lim T ^ob" 2 /-t 2 ] =92- &2k 2 > . 



(25) 
(26) 

(27a) 
(276) 
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For 82 > (] the wavefunction \& is connected with the Whittaker and the confluent series as 
follows 



or 



M/ = (C')- 1 / 2 M_^^ 1 (y^C) 



* = (C')- 1/2 (^C)^e-^/ 2 F(2 + J* ^) . 

2 4V02 



The series terminates if 



7 #1 



2 4^ 
The potential J7[C(x)] takes the form 



n, n = 0, 1, 2... 



2MC/[C(x)] 



g2C 2 +glC + ^0 + l 



+ 



(7 1 ^ 5A 



i? 2 ' 



where A = <t 2 — 4<72Co; here the function is given in the following implicit form 
±2(x - x ) = VR + y/coln\<Ti - 2t^\ + °"j_ ln| ^ + , 

2a/<72 t - y/(T2 

where i = ^p^ . 

It is supposed that i?(C) has no zeros within (0, 00), that is 



(28a) 



(286) 



(29) 



(30) 



(31) 



0"2 > 0, Co > 0, CTl > -2^/cq02 



(32) 



If two of the three parameters &%, 02 and Co turn to zero, Eq.(31) is easy to solve with respect to 
£ and the potential leads in this case to: a) a spherically symmetric oscillator well (or barrier) for 
o"2 = Co = pH ; b) a Morse potential for <7i = 02 = [12]; a Kratzer potential for <7i = cq = 
0,0]- 



We remark that in the work [15| an implicit potential occurs for which the Schroedinger 
equation is turned into a confluent HGM equation. This potential can be obtained from Eq.(30), 
if 02 = 0. The assumptions in |ll| may be considered equivalent to a transformation ((x) not 
depending on energy and therefore the results of |TI| can be easily obtained in our framework. 

Let us study now the definition interval of x. For ( —* 00 x goes to ±00 as i^y^C if °"2 / 0, 
and as ±V^lC, if 02 = 0, <7i 7^ 0, because in the latter case 



±2(x -x ) = Vr+-± + Vco^ \(t 2 \ 



(33) 



3 For S 2 = 0, * = (C')" 1/2 C 1/2 ^ao (V^C 172 ) (see the formula 2, 162 (la) in 0) 



(for c 2 = 0, o"i = x varies from — oo to +oo). For Q — > 0, x — > =F oo if Co 7^ 0, and x — > £0 if 
c = 0. 

For ci2 7^ 0, the zero of the energy is determined by the condition U(oo) = 0, i.e., 52 = 0. In 
this case, one gets for the energy spectrum the following equation 

01 



2n+ 1 



9i 



ho + 1- c kl . 



From this, it follows that 

a) for g\ > 0, ho + 1 > 0, there is no discrete spectrum, because 
0-1 



(34) 



2^2 



kr, 



ho + 1- c k% < 



I fen I < 



(35) 



b) for gi > 0, ho + 1 < 0, a finite number of discrete levels is possible if a\ < 0; 

c) for gi < the discrete spectrum has an infinite number of levels converging to zero. For 
o"2 = the discrete levels are determined by the equation 



l + ho-c kl= 2n + l + 9 \ 



(36) 

and for Co 7^ the discrete spectrum may have only a finite number of energy levels (the case 
o"2 = 0, Co = corresponds to an infinite number of levels). 

The wavefunctions of the discrete spectrum are given in terms of the Laguerre polynomials 
as follows 

Here, the normalization factor B n reads 

^ + ^ + (X + 2n + l)f) 
The integrals required for its calculation can be obtained from I v [O] : 

Iu = J e-^C-'lFi-n^-VhOM 
for v = 7 + 1,7. 

The potentials given by Eqs. (9) and (30) may have useful applications in the theory of 
diatomic molecules, in particular, for the calculations of important quantities such as the Franck- 
Condon factors and anharmonic constants. 



B,, 



.02 \ T(7 + n) 



rv. 



(37) 



(38) 



(39) 
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